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Abstract We report on the study of the momentum 
distribution of a one-dimensional Bose gas in an optical 
lattice. From the momentum distribution we extract the 
condensed fraction of the gas and thereby measure the 
depletion of the condensate and compare it with a the- 
orical estimate. We have measured the coherence length 
of the gas for systems with average occupation n > 1 
and n < 1 per lattice site. 



1 The one-dimensional Bose gas in an optical 
lattice 

A one-dimensional gas can be created in a trap when the 
confining potential restricts the motion of the particles 
to one dimension with the transverse motional degrees 
of freedom being frozen out. A cigar shaped harmonic 
trapping geometry is characterized by the frequencies 
uj± in the two strongly confining axes and ujax in the 
weakly confining axis. A kinematically one-dimensional 
situation is achieved when all particles occupy only the 
ground state in the radial directions, which implies that 
both the thermal energy ksT and the interaction energy 
/i have to be much smaller than the transverse energy 
level spacing. In general, one-dimensional systems ex- 
hibit increased quantum fluctuations of the phase, such 
that for a homogeneous ID gas Bose-Einstein condensa- 
tion is prevented. For trapped low-dimensional gases the 
cross-over to a Bose-Einstein condensate takes place at 
a finite temperature /cbTid = N'hujax/'^n{2N), where N 
is the number of particles in the ID system ^12]. The 
fluctuating phase alters the properties of the gas [3101|S1 
|B^ . One-dimensional quantum systems exhibit a wealth 
of fascinating phenomena whose explanations go beyond 
the mean- field description "5". 

One-dimensional trapped Bose-Einstein condensates 
were recently created and studied |^ using an optical 
lattice consisting of two mutually perpendicular stand- 
ing wave laser fields. In this geometry the optical lattice 



forms an array of one-dimensional tubes, each filled with 
a Bose-Einstein condensate. This experiment revealed 
the distinctively different excitation spectrum of a one- 
dimensional quantum system as compared to its three 
dimensional counterpart |1U| . In a previous experiment a 
Bose condensates was loaded into a two-dimensional op- 
tical lattices to study the coherence between the tubes. 
In that experiment the tunnel-coupling between adja- 
cent tubes was larger than the axial oscillation frequency, 
thereby an array of strongly coupled tubes was created 
|11|. Very recently the lifetime of one-dimensional gases 
created in an optical lattice were studied ^21- In elon- 
gated magnetic and optical traps a regime was accessed 
where a Bose condensate with fi < huj± coexisted with a 
three-dimensional thermal cloud |13II14| . In similar elon- 
gated traps studies of solitons ^21 and of enhanced phase 
fluctuations have been performed [TCItTT) . 

When the one-dimensional Bose gas is exposed to an 
additional optical lattice potential in axial direction the 
bosons may become localized in the minima of a periodic 
potential and the system can then be described by the 
Bose- Hubbard Hamiltonian |18II19| : 



H 






2 ^ 



(1) 



J denotes the hopping matrix element between neigh- 
boring lattice sites and determines the rate of which a 
particle disappears from lattice site i and tunnels to the 
adjacent lattice site j {hi and a\ are the annihilation and 
creation operators for an atom at lattice site i, respec- 
tively). The total tunnel coupling is J = z J with z being 
the coordination number of the lattice. The inhomogen- 
ity of the atom trap is characterized by e^ = ^{lujaxdY 
where m is the atomic mass and d the lattice spac- 
ing. The occupation number of lattice site i is denoted 
Tii and U is the onsite-interaction energy between two 
bosons on the same lattice site. The one-dimensional 
Bose-Hubbard Hamiltonian exhibits a transition from a 
superfluid phase to a Mott insulating phase for a ratio 
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{U/J)c — 2 OEIIHIEI, which we experimentally demon- 
strated recently |3- Due to stronger quantum fluctua- 
tions of the phase in a one-dimensional quantum sys- 
tem this value is lower than the corresponding critical 
value in three dimensions (U/J)c = 5.8 !18119i20|. Tun- 
ing the depth of the periodic potential changes the pa- 
rameter U/J, which leads to increased effects of inter- 
particle interactions and small filling of the lattice sites 
may turn the system into a gas of hardcore bosons j^fl 
I22| . In the inhomogeneous system Mott insulating re- 
gions with commensurate filling coexist with superfluid 
regions with incommensurate filling and the insulator is 
formed in a crossover type transition ^. 



2.2 Preparation of one- dimensional quantum gases 

Using the optical lattice we realize one-dimensional quan- 
tum gases ordered in an array of one-dimensional tubes. 
Two lattice axes are ramped to a fixed value V± = 
Vx = Vz = 30 Eji. In this configuration the transverse 
tunnelling rates Jx and J^ are small compared to the 
duration of the experiment. We therefore neglect tun- 
nelling between neighboring tubes and assume an array 
of isolated ID gases. An additional optical lattice is ap- 
plied along the symmetry axis with a potential depth 
Vax = Vy <^ V±. The potential depth Vax controls the 
value oi U/J. 



2 Experimental setup 



3 Results 



3.1 Momentum distribution 



2.1 Bose condensates in an optical lattice 



In the experiment, we collect up to 2 x 10^ ^^Rb atoms in 
a vapor cell magneto-optical trap. After polarization gra- 
dient cooling and optical pumping into the |F = 2, mp = 
2) hyperfine ground state the atoms are captured in a 
magnetic quadrupole trap. After magnetic transport of 
the trapped atoms over a distance of 40 cm the magnetic 
trapping potential is converted into the harmonic and 
elongated potential of a QUIC trap [22 • Subsequently, 
we perform radio frequency induced evaporation of the 
cloud over a period of 25 s. After evaporation we ob- 
serve almost pure three-dimensional Bose-Einstein con- 
densates of up to 1.5(0.2) X 10^ atoms. Following the 
condensation we adiabatically change the trapping ge- 
ometry to an approximately spherical symmetry with 
trapping frequencies of Wa; = 27rx 18 Hz, tjj, — 27rx20IIz, 
and a;^ = 27r X 22 Hz. This reduces the peak density by a 
factor of 4 and allows us to load the optical lattice more 
uniformly. 

The optical lattice is formed by three retro-reflected 
laser beams. Each beam is derived from a laser diode at 
a wavelength of A = 826 nm. At the position of the con- 
densate the beams are circularly focused to 1/e^-radii 
of 120 /im {x and y axes) and 105 /im (z). The three 
beams possess mutually orthogonal polarizations and 
their frequencies are offset with respect to each other 
by several ten MHz. We stabilize the lasers to a high- 
finesse Fabry-Perot cavity, thereby reducing their line 
width to ^ 10 kHz. In order to load the condensate into 
the ground state of the optical lattice, the intensities of 
the lasers are slowly increased to their final values us- 
ing an exponential ramp with a time constant of 25 ms 
and a duration of 100 ms. The resulting optical potential 
depths V^,j,,z are proportional to the laser intensities and 
are conveniently expressed in terms of the recoil energy 
^« = ^withfc=f. 



We have studied the momentum distribution of the gas 
by imaging the expanding atom cloud after being re- 
leased from the optical lattice. Prior to switching off 
the lattice the laser intensity of the axial lattice is in- 
creased to ~ 25Eii for 40 /iS. Then all laser beams are 
extinguished simultaneously within 5/j,s and the mag- 
netic trapping potential is switched off within 300 /iS. 
This procedure may be regarded as a projection of the 
wave function of the atoms onto the Bloch momentum 
states for a lattice depth of 25 E^. This method offers 
the advantage of controlling adiabatic processes during 
the switch-off of the lattice. These might affect the mo- 
mentum distribution of the released atoms '^ . 

For our one-dimensional gases the coUisional inter- 
action during expansion - which would give rise to a 
Thomas-Fermi type envelope of the density distribution 
in the y-direction - is small. This is related to the ex- 
tremely fast radial expansion of the clouds, which are 
initially confined to the harmonic oscillator ground state 
a± = y^h/{muj±) with a trapping frequency of lli± — 
27r X 33 kHz. The peak density decreases proportional to 
1 / y/l + (uj±T)^ which is much faster than the timescale 
for axial expansion that is given by the axial trapping 
frequency of ujax — 2t: x 80 Hz. This axial confinement 
is introduced by the gaussian waist of the laser beams 
forming the array of one-dimensional tubes. 

We have studied the resulting momentum distribu- 
tion after 15 ms of free expansion of the cloud. The atoms 
were imaged by resonant absorption imaging and we de- 
termine the optical density in the center of the cloud 
by averaging over a 62 /im wide section in the image. 
Fig. ^ shows the measured momentum distribution for 
a sample with an average filling of n — 1.2 atoms per 
lattice site. We calculate the average filling from the 
mean density of atoms which is derived known atom 
number and the measured trapping frequencies and in- 
cludes a modified coupling strength due to the localiza- 
tion of the atoms in the optical lattice potential ^25| . 
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momentum [tik] 

Fig. 1 Measured momentum distribution for A^ = 1.5 x 10^ 
atoms for different deptfi of the axial lattice potential. Top 
to bottom on the left hand site: Vax = 5, Vax = 7, Vax = 9, 
Vax = 12, Vax ~ 18 (in units of Eu). The inset shows an 
absorption image after 10 ms time of flight where the area of 
averaging is indicated. 



The atom number distribution at a given lattice site in 
the superfluid phase can be approximated by a Poisson 
distribution. Therefore, even with an average occupation 
number of n = 0.6 the probability of finding doubly or 
multiply occupied lattice sites is on the order of 10%. 



tight binding regime but for small depletion. Figure |21 
displays the measured coherent fraction together with 
the calculated values for the quantum depletion, which 
show good agreement in the range of validity of eq. ^ . 
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Fig. 2 Coherent fraction of the ID Bose gas as measured 
from the time-of-flight image. One-dimensional system with 
n — 1.2 (solid squares), one-dimensional system with n = 0.6 
(open squares), three-dimensional gas with n = 1 (triangles). 
Part of the experimental data are taken from |7|. The solid 
line shows the calculated quantum depletion according to 
equation ||5J using N^ = 40. 



3. 2 Quantum depletion of the condensate 

When the potential depth of the lattice is raised the 
strength of the atom-atom interaction increases. This 
is accompanied by quantum depletion of the conden- 
sate PD- The reduced condensate fraction of the sys- 
tem diminishes the contrast in the matter wave inter- 
ference pattern after the atoms have been released from 
the optical lattice. To extract the number of coherent 
atoms Ncoh from the interference pattern, the interfer- 
ence peaks at Ohk, ±2fik and ±Ahk are fitted by gaus- 
sians. Incoherent atoms, both due to localization of the 
atoms in the lattice or due to depletion of the conden- 
sate, give rise to a broad gaussian background centered 
at zero quasi momentum which dominates for higher Vax 
PS] . Taking this fit as a measure of the number of inco- 
herent atoms Nincoh, we calculate the coherent fraction 

fc = Tf — ^^# U\- The measured data are shown in 

figure 121 For the one-dimensional Bose gas in an optical 
lattice the quantum depletion is calculated by j26| 
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Here m is the atomic mass, m* in the effective mass in the 
lattice, uiD is the ID density, aiu is the one-dimensional 
scattering length 27 and N^ is the number of potential 
wells occupied by atoms. This expression holds in the 



3.3 Coherence length 

From the width of the central coherent momentum com- 
ponent the coherence length of the gas in the optical 
lattice can be inferred. A Bose-Einstein condensate ex- 
hibits a long coherence length which leads to a small 
width of the interference peak. Due to the inhomogene- 
ity of the trap the transition to the Mott insulator sets 
in at those positions, where the local density becomes 
commensurable with the spacing of the potential wells 
in the optical lattice. This process breaks up the con- 
densed cloud into smaller units and therefore reduces 
the coherence length of the system. In figure |21 we show 
the change of the coherence length as a function of the 
parameter U/J. For the one-dimensional gas with an 
average occupancy of n = 1.2 we observe a kink in the 
coherence length at [// J ~ 2 which indicates the onset 
of the Mott insulating phase. For the three-dimensional 
gas with fi — 1 we find the Mott insulator transition at 
U/J ~ 6, in agreement with the mean field theory. For 
very small one-dimensional systems with n = 0.6 we ob- 
serve that the coherence length of the sample starts to in- 
crease around U/J = 4. This is in qualitative agreement 
with the theoretical prediction that the Mott-insulator 
transition occurs for U/J ~ 2 only for a chemical poten- 
tial corresponding a density n = \ but at larger values 
for n < 1 0. 
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Fig. 3 Width of the central momentum peak as a measure 
of the phase coherence length of the system. One-dimensional 
system with n = 1.2 (solid squares), one-dimensional system 
with n = 0.6 (open squares), three-dimensional gas with n = 
1 (triangles). Part of the data are taken from ref. [7|. 
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4 Conclusion 

We have studied the momentum distribution of one- 
dimensional Bose gases in an optical lattice for different 
average occupation of the lattice sites. We find that the 
coherent fraction of atoms is independent of the atomic 
density when we tune the interaction strength U/J from 
the superfluid to the Mott-insulating regime. The change 
of the coherence length appear to indicate that the sys- 
tem enters the Mott-insulator for small atomic densities 
only for larger values of U / J . 
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Giamarchi, Corinna KoUath and Stefan Wessel and sup- 
port from SNF and QSIT. 
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